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Abstract. We present a detailed theoretical study of light scattering off superfluid 
trapped Fermi gas of atoms at zero temperature. We apply Nambu-Gorkov formalism 
of superconductivity to calculate the response function of superfluid gas due to 
stimulated light scattering taking into account the final state interactions. The 
polarization of light has been shown to play a significant role in response of Cooper- 
pairs in the presence of a magnetic field. Particularly important is a scheme of 
polarization-selective light scattering by either spin-component of the Cooper-pairs 
leading to the single-particle excitations of one spin-component only. These excitations 
have a threshold of 2 A where A is the superfluid gap energy. Furthermore, polarization- 
selective light scattering allows for unequal energy and momentum transfer to the 
two partner atoms of a Cooper-pair. In the regime of low energy (« 2 A) and low 
momentum (< 2A/(Tivf), Vf being the Fermi velocity) transfer, a small difference in 
momentum transfers to the two spin-components may be useful in exciting Bogoliubov- 
Anderson phonon mode. We present detailed results on the dynamic structure factor 
(DSF) deduced from the response function making use of generalized fluctuation- 
dissipation theorem. Model calculations using local density approximation for trapped 
superfluid Fermi gas shows that when the energy transfer is less than 2Ao, where Ao 
refers to the gap at the trap center, DSF as a function of energy transfer has reduced 
gradient compared to that of normal Fermi gas. 
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1. Introduction 

Cold atoms are of enormous research interest in current physics. The tremendous 
advancement in technology of cooling, trapping and manipulation pQ of atomic gases 
during 80 's and 90 's has enabled researchers to achieve a low temperature down to 
a few hundredth of a microKelvin. This led to the first realizations of Bose-Einstein 
condensation (BEC) in dilute gases of ultracold bosonic atoms about a decade ago. 
Predicted in 1924 by Einstein p. based quantum statistics of indistinguishable particles 
discovered by Bose [I], BEC in gaseous systems had long been thought a subject 
of mere academic pursuit beyond experimental reach because of the requirement of 
ultralow temperature which was unimaginable even two decades ago. The success in 
BEC is a breakthrough prompting researchers to look for experimental realizations of 
many other theoretical predictions of quantum physics using cold atoms. The most 
remarkable property of such atoms is the tunability of the atom-atom interaction over 
a wide range by an externally applied magnetic field or other means. This provides an 
unique opportunity to explore physics of interacting many-particle systems in a new 
parameter regime. In this context, the focus of attention has been now shifted to 
cold atoms obeying Fermi-Dirac statistics. Since fermions are the basic constituents of 
matter, research with Fermi atoms under controlled physical conditions has important 
implications in the entire spectrum of physical and chemical sciences. In particular, it 
has significant relevance in the field of superconductivity jSlEj- 

The quantum degeneracy in an atomic Fermi gas was first realized by Jin's group 
[7j in 1999. Since then, cold Fermi atoms have been in focus of research interest in 
physics today. In a series of experiments, several groups jHJ 13 EE EH 112 UH] have 
demonstrated many new aspects of degenerate atomic Fermi gases. In a remarkable 
recent experiment, Ketterle's group has realized quantized vortices as a signature of 
Fermi superfluidity in a trapped atomic gas. Two groups-Innsbruck ^H] and JILA 
[To] have independently reported the measurement of pairing gap in Fermi atoms. 
Furthermore, Duke and Innsbruck groups ^3 have measured collective oscillations 
which indicate the occurrence of superfluidity [JJJ]. One of the key issues in this field 
is the crossover [2H1 EU E2] between BCS state of atoms and BEC of molecules formed 
from Fermi atoms. Several groups have demonstrated BEC [2S] of molecules formed from 
degenerate Fermi gas. There have been many other experimental [21] and theoretical 
investigations [2H] revealing many intriguing aspects of interacting Fermi atoms. 

The analysis of response of Cooper-paired Fermi atoms due to external perturbation 
(such as photon or rf field) is important for understanding the nature of atomic Fermi 
superfluid. A method has been suggested to use resonant light j2Ej to excite one of 
the spin components into an excited electronic state and thereby making an interface 
between normal and superfluid atoms. This is analogous to superconductive tunnelling 
which has a threshold equal to the gap energy A. This has been recently implemented 
(albeit using rf field) [THJ [22j to estimate gap energy. There have been several other 
proposals [2U [2E| for probing pairing gap. 
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Our purpose here is to calculate response function of superfluid Fermi gas due to 
stimulated light scattering that does not cause any electronic excitation in the atoms. 
We particularly emphasize the role of light polarization in single-particle excitations 
which have a threshold 2 A. We present a scheme by which it is possible to have 
single-particle excitation in only one partner atom (of a particular hyperfme spin 
state) of a Cooper-pair using proper light polarizations in the presence of a magnetic 
field. This may lead to better precision in spin-selective time-of-flight detection of 
scattered atoms. Furthermore, spin-selective light scattering allows for unequal energy 
and momentum transfer into the two partner atoms of a Cooper-pair. This may be 
useful in exciting Bogoliubov-Anderson (BA) phonon mode of symmetry breaking by 
making small difference in momentum transfers received by the two partner atoms from 
the photon fields. A number of authors [SH EDI E] have theoretically investigated 
Bogoliubov-Anderson (BA) mode pffi ESI E3| in fermionic atoms as a signature of 
superfluidity. BA mode constitutes a distinctive feature of superfluidity in neutral Fermi 
systems since it is associated with long wave Cooper-pair density fluctuations. However, 
experimental detection of this mode is a challenging problem. 

We present a detailed theoretical analysis of the response function of Cooper-paired 
atoms at zero temperature due to light scattering. The stimulated light scattering we 
discuss here is similar to Bragg spectroscopy used by Ketterle's group for measuring 
structure factor of an atomic BEC The response function we derive is applicable 
for most general case of polarization-selective single-particle excitations for unequal 
(or equal) momentum as well as energy transfers to the two partner atoms of a 
Cooper-pair. We develop the theoretical framework for stimulated light scattering 
off Cooper-paired Fermi atoms following the method used for describing Raman 
scattering in superconductors jSHl HH] • We use standard Nambu-Gorkov formalism of 
superconductivity [SHI HH] to calculate the response function taking into account the 
vertex correction due to final state interactions. We deduce dynamic structure factor 
(DSF) from the response function applying generalized fluctuation-dissipation theorem. 
We present detailed analytical and numerical results of our calculation of DSF of trapped 
superfluid Fermi gas of atoms using local density approximation. The inhomogeneity of 
trapped gas has a role in distinguishing the DSF of superfluid gas from that of normal 
gas. When the energy transfer is smaller than 2Aq where Ao is the gap at the trap 
center, the DSF of superfluid g function of energy transfer shows much reduced 

gradient in comparison to that of normal gas. This is because of the fact that the gap 
A has an inhomogeneous distribution gradually vanishing at the edge of the trap. 

The paper is organized in the following way. In the following two sections, we define 
bare vertex in light scattering and response function, respectively. In the fourth section, 
we discuss stimulated light scattering in two-component 6 Li Fermi atoms in the presence 
of a magnetic field. We next describe in detail the method of vertex correction in light 
scattering off Cooper-paired Fermi atoms. In the sixth section, we discuss our analytical 
results followed by description on numerical results in the seventh section, and then we 
conclude. 
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2. Bare vertex in light scattering 

To begin with, let us consider an elementary process of photon scattering by a neutral 
atom. Let the atom's initial and scattered electronic state be denoted by A and B, 
respectively. The frequencies of the incident and scattered photon are represented by 
uj\ and u>2, respectively. According to second order perturbation theory, the strength of 
scattering is given by Kramers-Heisenberg formula 



iBA = 8ABh ■ £2 T 



(P-e 2 )MP-ei)/A , {P-£i)Bi{P-e2)iA 



1) 



LO IA - LO X UJ IA + UJ 2 

where / denotes all the intermediate atomic states that can be coupled to the initial 
and final atomic states A and B by the incident and scattered photon fields. Here 
p and m e are the momentum and mass of the valence electron of atom, ei( 2 ) denotes 
the polarization state of the incident (scattered) photon, ojja is the atomic frequency 
between the states / and A. The atomic transition (A — > B) probability and the 
differential scattering cross section of photons is proportional to |7| 2 [4*0] . It should be 
mentioned that 7 does not depend on the momentum transfer q associated with the 
scattering, but it is sensitive to light polarization directions. Let us now consider the 
particular case: A = B that is, before and after the scattering, the atom remains in the 
same electronic state. Then, making use of the completeness of the intermediate states 
/, one can rewrite the term 61.62 as jlO] 

ei-e 2 = ~~r — [(p-^)ai{pA)ia + (pA)ai{p^2)ia] , (2) 
m e n j io 1 a 

Further, let us assume u>i ~ uj 2 — ujja, that is, the incident as well as scattered light 
fields are in near resonance with the atomic frequency. In such a case, the second term 
within the third bracket on the right hand side (RHS) of Eq. (0) is much smaller than 
the first term, because energy denominator of the second term is of the order of optical 
frequency while that of the first term can be chosen to be smaller by several orders of 
magnitude. Thus, neglecting the second term, the bare vertex can be written as 



1 ^o;i(p.e 2 ) Ar (p.ei) /A 
Iaa — 7 \— • 3 

Next, using electric-dipole approximation and the fact uji/ujja — 1, one can express 

n -iv^(di/.4)(dA/^i) 
j h [ujai - oj{) 

where d^/ is the transition dipole moment between the states A and /, £ j = S^i is the 
electric field and 

e 2 £ 1 8 2 

^0 = = (5) 



with ni being the number of incident photons. 
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3. The response function 

To define response function of fermionic atoms due to an applied laser field, we use the 
second-quantized operator (V^aJ.jJ which describes the annihilation(creation) of an 
atom with hyperfine spin a and center-of-mass momentum k. These operators satisfy 
fermionic algebra. The effective atom- field hamiltonian is H e ff = Hq + Hi, where 

hk 2 

H o = J2 H u k ~ 6)c4 k a CT)k , uj k = — , (6) 

<7,k ~~ 

with 5 = u)\ — 0J2 being the frequency- difference between incident and scattered photons. 
We assume that, except the center-of-mass momentum, the spin or any other internal 
degrees of atom does not change due to light scattering. By treating light fields 
classically, the effective interaction hamiltonian can then be written as 

Ei = hQ 7<xo-«i,k+q a <T,k + H.c. (7) 
<t, k 

where q is the momentum transferred to the atom due to photon scattering and 7 CTCT 
represents the bare vertex corresponding to the ground hyperfine spin magnetic quantum 
number a. 

Now, one can define the density operators by = Z)<j,k a t,k+q a °-,k and 

P ( q ] = 7<T<Tai,k+qOff,k (8) 
k,cr 

One can identify the operator p^ as the Fourier transform of the density operator in real 
space. The scattering probability of incident particles (photons in the present context) 
is related to the response or susceptibility 

X(V,T-T') = -(T T [pW(T) P { 2l(T'))) (9) 

of the target system by which the incident particles are scattered. Here (■ ■ •) means 
thermal averaging and T T is the complex time r ordering operator. The Fourier 
transform of this susceptibility is 

1 f T 



X(q,w n ) = r/ dre^ T X (ci,r) (10) 

Z J-T 

where T is the temperature and u n = lisnT is the Matsubara frequency with n being an 
integer. The scattering cross section is proportional to the generalized dynamic structure 
factor which can be obtained by the generalized fluctuation-dissipation theorem through 
the analytic continuation of x(q, w„) as 

S(q,u) = --[l+n B M]Im[x(q,z = u; + i<5)]. (11) 

7T 

We define the following polarization matrix element: 

U tJ ( q ,r-r') = -(T T [p?(r)p^(r')]) (12) 

where i,j = 7,0. The polarization bubble IT 77 is nothing but the susceptibility 
x(q, r — r') of Eq. ©. The dynamic structure factor is thus related to this polarization 
term by fluctuation-dissipation relation as expressed in Eq. The spectrum of 

density fluctuation is proportional to the dynamic structure factor which can also be 
defined as the Fourier transform of the two-time density-density correlation function. 
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Figure 1. A schematic level diagram for polarization-selective light scattering in 
two-component Fermi gas of 6 Li atoms 



4. stimulated light scattering in two-component Fermi atoms 

We would like to study stimulated light scattering in two-component Fermi atoms. In 
particular, we consider trapped 6 Li Fermi atoms in their two lowest hyperfine spin states 
| g) 1 =| 2S 1/2 ,F = 1/2, mp. = 1/2) =|t) and | g) 2 =\ 2S 1/2 , F = l/2,m F = -1/2) = \i). 
For simplicity, the number of atoms in each spin component is assumed to be the same. 
However, a mismatch in number densities of the two spin components may lead to 
interior gap superfluidity |4*T) 14*2*] in a Fermi gas of atoms. An applied magnetic field 
tuned near the Feshbach resonance (~ 850 Gauss) results in splitting between the two 
spin states by ~ 75 MHz 43J, while the corresponding splitting between the excited 
states | e)i =| 2P 3/2 ,F = 3/2, m F = -1/2) and | e) 2 =| 2P 3/2 ,F = 3/2, m F = -3/2) is 
~ 994 MHz [Ej. 

Figure 1 shows the schematic level diagram for stimulated light scattering by two- 
component 6 Li atoms. Two off-resonant laser beams with a small frequency difference 
are impinged on atoms, the scattering of one laser photon is stimulated by the other 
photon. In this process, one laser photon is annihilated and reappeared as a scattered 
photon propagating along the other laser beam. The magnitude of momentum transfer 
is q ~ 2ki sin(#/2), where 9 is the angle between the two beams and is the momentum 
of a laser photon. Let both the laser beams be cr__ polarized and tuned near the transition 
I 9)2 — H e )2- Then the transition between the states | g)\ and | e) 2 would be forbidden 
while the transition | g)i — >\ e)\ will be suppressed due to the large detuning ~ 900 
MHz. This leads to a situation where the scattered atoms remain in the same initial 
internal state | g) 2 . Similarly, atoms in state | g)i only suffer scattering when two 
a + polarized lasers are tuned near the transition | g)i —>-\ 2P 3 / 2 ,F = 3/2, m F = 3/2). 
Thus, it is possible to scatter atoms selectively of either spin components using circularly 
polarized lasers in the presence of magnetic field. Under such conditions, considering a 
uniform gas of atoms, the effective laser-atom interaction Hamiltonian in electric-dipole 
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approximation can be written as 

Hi = tiQo 7<r<T<4,k+q a o-,k + H.c. (13) 

k,<x=U 

If a refers to \[) then 

O-i (d22-g 2 )(d 22 .gi) 

i=l,2 « (^22 - Wi) 

where da is the transition dipole matrix element between the ground | g)i and the excited 
| e)i states. Similarly, if a is ||) then the subscript "22" should be replaced by "11". 
In writing the above vertex term, we have also assumed that both the laser beams are 
of almost equal intensity. For both the laser beams having <r_ polarization tuned near 
I 9)2 — H e} 2 as in Fig. 1, one finds 7^ » 7^. On the other hand, in the absence 
of magnetic field (or in the presence of a weak magnetic field), the hyperfme magnetic 
sub- levels of the ground and excited state would be degenerate (or nearly degenerate). 
In such a case, irrespective of whether both the laser beams are unpolarized or equally 
polarized, we have 7^ ~ -yn- 

5. light scattering in Cooper-paired Fermi atoms: Vertex correction 

To study light scattering in Cooper-paired Fermi atoms, we apply Nambu-Gorkov 
formalism that uses the four Pauli matrices 

to= {1 ?)' n= {°i I) r2= (! 7)' r3= (i -1) (15) 

The vertex equation is |H] 

F(k + , + 0^r 3 G{k' + )T{k' + , k'_)G(k'_)r 3 V(k, k'), (16) 

where k± = k±q/2 and k = (k, ko) is the energy-momentum 4- vector whose components 
are k 3 = £/. and k^ = ik®. In pairing approximation, the Green function can be expressed 
in a matrix form as 

G(k) = W^WTTs ■ < 17 > 



where Ey. = y ^ + and = — ji with ek = h k 2 /{2m). The bare vertex 

'"(7 4)- 

Using Pauli matrices To and t 3 , this can be rewritten as 

7 = 7o^o + 73^3, (19) 
where 70 = [711 — 7u]/2 and 73 = [711 + 7j.jJ/2. The susceptibility is given by 

XM = J ^Tr[7 k G(k + )r(k + ,k_)G(k_)] (20) 
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Figure 2. The upper part shows two irreducible four- vertex diagrams of stimulated 
light scattering off noninteracting normal Fermi gas of atoms. The thin dashed lines 
with arrows represent the incident laser field (with electric filed Ei) and the thick 
dashed lines refer to the emitted photon stimulated by another laser field (E s ). The 
operators d.E^ and d.E s act at the vertex pairs (1,4) and (2,3), respectively; where 
d is the transition dipole moment between the ground | g\) (| (72)) and the excited 
I e i) ( I e 2) ) state. The vertex pair (1,2) can be replaced by an effective single vertex 
where the operator 7 CTCT p q acts, where 70-0- is given as in Eq. (|14fl . Similarly, the pair 
(3,4) can be combined to form an effective vertex. Thus, the four vertex diagrams 
effectively reduce to bubble diagrams as shown in the lower part. Note that the role 
of incident and scattered fields can be reversed, since an atom can absorb a photon 
from the laser mode marked "E s " and emit into the mode marked "Ej". By treating 
laser fields classically, the effective vertex operators can be expressed only in terms of 
atomic Fermi operators as in Eq. i|13|) . 




Figure 3. An irreducible four- vertex diagrams of stimulated light scattering in an 
atomic Fermi superfiuid when the quasi-particles are assumed to be noninteracting. 
The double lines with arrows represent Nambu propagator for Cooper pairs. As in 
Fig 2, the four vertex diagram can be effectively represented by a two-vertex bubble 
diagram. 




Figure 4. Diagrammatic representation of vertex equation 
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5.1. vertex equation and its solution 

To solve the vertex equation, let us expand the vertex function in terms of Pauli matrices 
as 

T(k + ,k^=j2T^(k, q ,u)T l . (21) 

4 = 

Using Eqs. (fT7|) and (J2*T|) in Eq. (fTKj) . we can write 

r(k, q, uj) — % + i J ^_V(k, k') {k/f + A ^ (fc , 2 - ^ g ^ (22) 
where A; 2 = fe| + k\ = £ k — k% and 

yo = W +0 e k >_ + k'_ Q e k , + ]rw _ iA k ,\e k , + - e k ,_]T^ + A k/ [k' +0 + fc'_ ]r« 

+ [e k/+ e kL +k' +0 ki + A 2 k/ ]T^, 

Vl= - i[k' +0 e kL - k'_ e K }T^ - A k ,[e K + e kL ]T^ - A ¥ [k' +0 + k'_ ]T^ 
+ [e K e kL -k' +0 k'_ -A 2 ]T^\ 

y 2 = - iA k ,[k' +0 - k'_ ]T® - iA k ,[e k , + - e w _\Y<® - i[e k , + k'_ - k' +0 e kL ]T^ 
+ {e K e kL - k' +0 k'_ + A 2 k ,p 2 \ 

y 3 = A k ,[e K + e kL ]rW - iA k ,[k' +0 - ^ ]r (2) + [e y+ e kL + k' +o k'_ - A 2 k ,]T^ 
+ [e kf+ k'_ + k' +0 e kl _}V^. 

In writing the above equations, we have assumed A k± ~ A k . Further, we can write 
e k± = £ k ± v fc .p 9 /2 + e q , where v fe = hk/m, p q = hq and e q = p 2 q /{2m) 

Before performing the integration of Eq. ([22)1. we note that the dominant 
contribution to the integral comes from fc-values near ~ 0, that is, h 2 k 2 /(2m) ~ fx. 
Hence we can approximate 

r d 4 k r d 3 k r r dk 3 dk 

J J wy 6 ^J J— ^ 

where k 3 = denotes the third component of energy-momentum 4- vector k = (k, ko). If 
the potential V(k, k') is separable in two variables k and k', then Eq. (JT5J1 is analytically 
solvable. Let us, for simplicity, replace V(k, k') by the well-known mean field potential 
V m f = ga s (where g = Anh 2 / '(2m)) which is expressed in terms of s-wave scattering 
length a s . By doing so, we are basically considering the weak-coupling case. However, 
within mean-field approximation the strong-coupling limit may be accessed by first 
renormalizing the BCS mean-filed interaction and then taking the limit a s — > ±oo as 
will be discussed later. 

With the assumption of a ^-independent gap A, the double integrations on k and 
&3 then resemble to those appearing in relativistic equations in QED and so can be 
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carried out analytically by Feynman's method The angular integration is left to 
the last. There are basically two types of integrals: 

A 2 dk dk 3 



HQ) 



(h\ + A 2 )(P + A 2 ) 



hJ 



3,4 



(24) 



(25) 



(k\ + A 2 ) (A; 2 + A 2 ) 

These integrals are explicitly calculated in Ref. [46] using Feynman's method of 
parametrization. For completeness, we here reproduce the method of calculation. The 
terms which are odd in k will not contribute to the integration and so those terms can 
be omitted. Substituting k = k — (q/2 — qx) where x is a parameter varying between 
to 1, the integral of Eq. (J2"4]) can be reexpressed as 



dx 



A 2 dkr)dki 



• (26) 

[k 2 + A 2 + q 2 (x - x 2 )] 2 

The fccr m tegration can be carried out by residue method of complex integration. The 

A 2 + q 2i 



pole is k 



k\ + L, where L 



(x 



x- 



Since L has infinitesimally 



negative imaginary part, the pole lies in the lower half of the real axis. The residue 
is — [4(^3 + L) 3 / 2 ] -1 . After performing k 3 - and x-integration , one obtains the result 
I(q) = /(?)/2, where 

arcsin j3 



6 2 



10 



v v fc-P 9 



(27) 



(2f 



4A 2 

The /^-integration in Eq. (}26|) is divergent, therefore a cut-off frequency uj c is required 
as the upper limit of integration. After having performed the integration, the vertex 
terms rw can be expressed as 

r d 3 ~k 

T^= l0 (k)-ga s — —Sfa) 



r (i) 



^(vfc-Pg)(l - f) 

U 2 - (Vfc.Pg) 2 
d 3 k 



p(3) 



,- Vfc-Pg/ 

2A 



F (2) + 



-ga. 



= - 



ga 



(2tt) 3 
d 3 k 



(2tt) 3 

Vfc-Pg/ 

2A 



ln iAT + (/?2 



1)/ 



K-p g ) 2 (i - f) 

U 2 - (Vfc.Pg) 2 
r (l) 



(29) 
(30) 



r (o) + {ln ^ + ^ /}r (2)_.^ r ,,, 



IAI 



2 A 



(31) 



r (3) 



d 3 k 



73 (A;) - ga, 

UJ f T (2) + (Vfc-Pg) 2 -^V r ( 3 ) + ^(Vfc-Pg)(l - f) 



—l- 



2A 



or 



(Vfc.Pq 



r (o) 



IVfc.Pg 



(32) 
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Since is decoupled from all other vertex terms including the bare ones (7$), we 
can set = 0. Using the expansion of Eq. ([21)1. the susceptibility can be written as 



2(r(°)- 7 o) 



X 



ga 

"\2 



7o + 2 



d 3 k 



5(6 



(Vfc.Pg) 2 - U 2 f 



p(3) 



,(0) 



73- (33) 



We note that the dressed part of is proportional to the momentum transfer q, 
therefore we have ~ 70 (k) in the low momentum transfer regime, that is, for 
q « where £ = hvp/(2A) is the BCS coherence length. Introducing the variable 
z = cos 9, where 6 is the angle between and p q , we can drop all the terms odd in 
z in the above equations, since upon integration over z those terms vanish. Thus 1^°) 
also becomes decoupled while r( 2 ) and form only two coupled equations which can 
be analytically solved. 



5.2. gap equation 

The gap equation can be obtained from Eq. (|3"T]) by setting q and uj equal to zero and 
replacing by the gap parameter A. The resulting equation reads 



A = -ga s J 



d 3 k r U uj c 



(34) 



The cut-off frequency u c has been introduced ad-hoc to tackle the divergence problem 
for the time being. This needs to be eliminated by the method of regularization. To this 
end, we here recall that in carrying out the various momentum integration, we made 
an approximation: the integration was restricted near the chemical potential (which is 
nearly equal to Fermi energy in the weak coupling regime). To restore the actual gap 
equation, we here remove this approximation and let u c — > 00 and thus obtain 

1 1 r d 3 k 1 

"2 



ga s 



(35) 



The gap defined by this equation is however, divergent. To remove this divergence, we 
define regularized mean-field coupling by subtracting from the right hand side of Eq. 
fl33|) the zero field contribution (i.e., A = and fi = 0). The resulting gap equation is 

d 3 k I I 



1 

ga s 



(27T) 



a + a 2 



(36) 



which yields convergent results. In the weak-coupling regime (|a s |/cp « 1), // ~ oc 
n 2 / 3 . The strong-coupling regime (|a s |fci? > 1) may be accessed by simultaneously 
solving for the interacting chemical potential fi from the single-spin BCS number- density 
equation 
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Figure 5. Gap A (in unit of ep) is plotted as a function of the dimensionless mean- 
field interaction parameter kp \a s \ on semi- logarithmic scale. The inset shows the same 
plot for small interaction parameter on linear scale. In the limit fci?|a s | — > oo, the gap 
saturates at a value 0.68 ep. 



This approach of solving the regularized gap plus the number equation to access strong- 
coupling regime within the simple mean-field framework fails to account for pairing 
fluctuation effects which are particularly significant near T c in the strong-coupling 
regime. However, far below T c , the correction due to the pairing fluctuation is very 
small as shown in Ref. [21]. The two coupled Eqs. (|36j) and pTj) admit analytical 
solutions which are obtained by Marini et al. jUj for the entire range of the parameter 
a s kp starting from weak interaction (a s kp — > ±0) to the unitarity limit (a s kp — > ±oo). 
In the unitarity limit, the solutions provide fj, = 0.59ep and A ~ 1.16//. For convenience 
in solving the two coupled equations numerically, we rewrite the equations in terms of 
the two dimensionless scaled variables x = k/k^ and y = A//i as 



2vr 



X 



x- 



dx 



(38) 




x 



X 



1)2 +y 2 



dx 



(39) 



where k^ = y/2mfi/h. We have set a s = — \a s \. Calling the right hand side of Eqs. (J3*%|) 
and (|3*9*j) as I\ and I2, respectively; eliminating k^ from both the equations, we obtain 
2tt h 

kM ~ W 1 *' [ } 

For given values of the parameters kp and \a s \, the Eq. ()4())1 can be solved for y. Then 
substituting this solution into Eq. (|3*8*j) . one evaluates \i and so also the gap A = \iy 
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Figure 6. Chemical potential \i (in unit of ep) is plotted against parameter /sp|a s |. 
In the limit fc_p|a s | — > oo, \x saturates at a value 0.59 ep. In the limit fci?|a s | — > 0, fi 
goes to unity. 



5.3. solutions 

Now, to write down the solutions of the various vertex terms and the susceptibility 
X is straightforward. Let k s = N(0)ga s , where iV(0) = (ir 2 ?i 2 )~ 1 mkF represents the 
single particle density of states near the chemical potential. The various vertex terms 
can be expressed as 

r <3> = (4i) 

p(2) _ iv(f) -p(3) (A2) 

r(0) = TT^y <«> 
.i { T"! 2 </• (45) 

^ - (Vfc.Pg) 2 

B= (v tPa f(l-/) 

The symbol (X) implies average of a function X over the chemical potential surface: 
(X) = [X(O)]^ 1 / d 3 k5(e fc )X, since X is an even function of z — cos9, we have 



and 



Here 



and 
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(X) = (1/2) Jg 1 X(kp, z)dz. Making use of these vertex terms, the susceptibility can 
be written as 

ll{B) 



X (q,u;)=2iV(0): 



2N(0) 



' 2 K s -f 3 F 2 

jZp 

' 3 1 + K S F 



(47) 



1 + K.(B) 

We drop the second term inside the third bracket which leads to small corrections due to 
Landau-liquid-like behavior without adding any significant qualitative effect. Further, 
for k s (B) « 1, we have 



X (q, u) = 2N(0) [ll (B) + 7 2 fJ (48) 
6. dynamic structure factor 

The dynamic structure factor is obtained from the response function x y i a analytic 
continuation of energy uj —>■ uj + i0 + . By means of generalized fluctuation-dissipation 
theorem as embodied in Eq. (jllj) . in the zero temperature limit the dynamic structure 
factor is related to the imaginary part of the density response function x y i a analytic 
continuation of energy u — > uj + i0 + as 

5(q, u) = --Im[x(q, u - uj + i0 + )].. (49) 

7T 

The key function here is /(/?) of Eq. (|2*7|). where (3 is given in Eq. (|2*5j) . As u; — > co> + i0 + , 
(3 — > (3 + i0 + . We have the following analytic properties of f{/3): 

m = m + p > 1 (50) 

where 

arcsinliA//? 2 — 1 



- evs^r (51) 

The use of Eqs. (jUJ, (BHJ) and (@Sl in Eq. (gBJ) which, along with Eq. flHDJ), on 
being substituted in Eq. (|4*9*j) leads to the result 

%f) =7o^ (q,^) + 7f£ 3 (q,o/) (52) 



where 



S 3 (q,uj) = N(0)^-/ M * \ - N(0) ^ ' 



4A 2 \/5 3 VF^T/ v 7 4A 2 |(/3 2 /)| 2 
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7. Analytical results and discussions 

Equation (|H2j) gives an expression for dynamic structure factor of a homogeneous Fermi 
superfluid when the excitations are of single-particle type for the parameters satisfying 
(3 > 1. Different amount of energy transfers (or excitations) to the two constituent 
partners of a broken Cooper-pair can be made by appropriately selecting the polarization 
states of the exciting two laser beams and tuning their frequency from the excited atomic 
state in the presence of a magnetic field. This fact is taken into account in the expression 
of (|52|h because any nonzero value of the term 7 means unequal excitation of the two 
partners. For instance, two extreme cases can be mentioned: Case-I: For unpolarized 
light in the absence of magnetic field, equal amount of energy transfer occurs to the two 
partners resulting in 70 = 0; Case-II: On the other hand, for circularly polarized light 
in the presence of strong magnetic field, we have 7q ~ 7! meaning only either partner 
can be excited. We will present our numerical results for these two extreme cases. To 
compare our results with the known results for normal Fermi system in the limit A — > 0, 
we will use in Case-I the limit 7^ ~ 7^ — > 1 meaning 70 — > and 73 — > 1. In Case-II, 
we will use the limit 7^ ~ jh — > 1 implying that 73 = —70 — > 1/2. Intuitively, 
one may understand that the Case-II would be significantly different from Case-I both 
qualitatively and quantitatively. In the Case-II, upon receiving an energy u (> 2A) 
from an incident photon, one partner of a Cooper-pair moves out of the Fermi sphere, 
while the other partner remains within the Fermi sphere. Let us consider an elementary 
process of single photon scattering by a Cooper-pair. Suppose, the Cooper-pair consists 
of an atom A having spin j and momentum k and another atom B with spin j and 
momentum — k. When this Cooper-pair is broken due to stimulated scattering of cr_ 
polarized photon in a situation like Case-II, atom A will move out of the Fermi surface 
as an excited quasi-particle with momentum k + q with certain probability given by 
BCS correlation and atom B will have certain probability of remaining within the Fermi 
sphere moving as a quasi-particle with momentum — k. Thus, only one partner of the 
Cooper-pair will contribute to the intensity of scattered atoms reducing the strength of 
the density fluctuation spectrum compared to that of Case-I. However, there could be 
some advantage in detecting the scattered atoms in Case-II by spin-selective time-of- 
flight measurement technique as we will discuss later in the concluding section. 

7.1. Case-I: Leading approximations 

In this case, we have 70 = 0. In the limit 73 — > 1, 



which is given by Eq. (JS3J). For (3 > 1, in the leading approximation in terms of (3 
this reduces to the form 




(55) 




(56) 
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which is devoid of any vertex correction. The same expression can be derived by taking 
p(3) _> -y 3? p(o) _> ry anc [ p(2) _^ q meaning that we use bare vertex only This is 
also obtainable from the static BCS- Bogoliubov mean-filed treatment as shown in the 
appendix. Because of the absence of vertex correction, it violates the Ward identities 
[IB] that guarantee the conservation of total particle number and the obeyance of the 
continuity equation. 

To perform the integration over z in Eq. (|56p. it is convenient to change the variable 

into 

x = , P « VkZ (57) 
Vuj 2 - 4A 2 V ; 

The condition (3 > 1 implies x < 1. Then the Eq. (JHH)) can be expressed as 
M , , .x 2iV(0)A 2 r-o 



5f a (q,^) = ^ / 



WPg^F JO (1 - JX 2 ) 3 / 2 V / T^X 2 " 

(5 

where j — 1 — 4A 2 /cj 2 and 



xq = Min 



y/u 2 - 4A 2 



(59) 



For 2A < lu < J (p q v F ) 2 + 4A 2 , we have x = 1 and the result is 



Sr( q ,u J ) = ^E(j), (60) 

where E(j) is the complete elliptic integral. Note that in the limit A — > 5} ead (q, a;) 
reduces to the form N(Q)u / (2p q vp) which is same as that of a normal quantum fluid of 
noninteracting quasi-particles within the energy range < ui < Vpp q jlHj- The dynamic 



structure factor reaches a maximum at loq = J (p g VF) 2 + 4A 2 . As lo increases above u>o, 
Xq decreases below unity and hence the integral in Eq. (|58j) decreases. 

In view of the forgoing analysis, we now verify how far f-sum rule is fulfilled by the 
dynamic structure factor as given by Eq. (|56p. To this end, we separate the integral 
over energy in the sum rule 

u>S(q, uS)duj — / ■■■doj+ •■•did. (61) 

Since Eq. (I56|) holds good for uo > 2 A, the first integral appearing on the right hand 
side of Eq. (JBTj) results in 

r"o N(0) r"o 



h= I ~ uS(q,u)dw = ipZL I uj 2 E{j)duo (62) 

J2A 2p q VF J2A 

where we have used the Eq. ([60)1 . In the limit A — > 0, E(j) — > 1 and so we obtain 

= N(0)( Pq v F ) 2 = ^pl_ = Nfi 
1 6 3yr 2 2m 2m 1 ' 

where we have used N(0) = {-Kh)^ 2 mkp. Here N represents the total number of particle 
per unit volume. The second integral on the RHS of Eq. (JoTj) is much smaller than 
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the first one. Thus, we find that in the limit A — > 0, or alternatively, for p q vp » 2 A 
and uj » 2A, that is, for large momentum and energy transfer, the dynamic structure 
factor as given by Eq. (j56|l approximately satisfies the f-sumrule. In this context, it may 
be worthwhile to mention here that for evaluating gap energy from the measurements of 
the scattering cross section of the light-scattered atoms released from a trap, large 
momentum transfer is indeed required to distinguish the scattered atoms from the 
un-scattered ones I50j . For single-particle excitation (/3 > 1) with small energy 
transfer, this leading approximation is not valid and the second term on the RHS 
of Eq. (p>fj) makes significant contribution resulting from vertex correction. We will 
show in the appendix that the DSF in leading order approximation is obtainable from 
BCS-Bogoliubov mean-field treatment that does not take into account final state (quasi- 
particles) interaction. 

7.2. Case-II: Leading approximations 

In this case 7^ — > and 7^ — > 1 implying 7^ ~ 73 7^ 0. Let us use 7q ~ 73 — > 1/4. 
Then we have 

Su{q,u) = . (64) 

The angular integration in Eq. (p)3*j) can be conveniently performed using the x-variable 
as already introduced in Eq. (joT|) . Explicitly, this takes the form 

2iV(0)A 2 p jx 2 
S (q,u}) = / dx- ———== 65 

OJPqVp JO (1 - JX 2 ) 3 / 2 Vl - X 2 



where j and xq are already defined above. For 2A < lu < y (p g VF) 2 + 4A 2 , xo = 1 and 
the result is 

S (q,u) = 2iV(0)A2 ^ 2 F 1 (3/2,3/2;2,j) (66) 
up q v F 4 

2Fi(a, b; c, d) is the hypergeometric function. In the limit A — > 0, 7rJ ^~ J ' > 2 -Pi ~~ ^ 1> hence 
we obtain 

«->-^ 

which again coincides with the form of the DSF of normal fluid within the specified 
parameter regime. 

In passing, we reemphasize that the leading order approximations are valid for 
j3 » 1, that is, for large energy transfer. In this limit, j — > 1 and so DSF's in both the 
cases tend to become equal. All he results followed from leading order approximations 
can also be obtained with BCS-Bogoliubov mean-field approximation without any vertex 
correction as illustrated in the appendix. 
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7.3. Bogoliubov- Anderson mode 

Now let us consider the case < (3 « 1, that is v^.pg < u> « 2A. In this case, the 
second term in Eq. (|4"4"|) dominates over all other terms. This term leads to Bogolibov- 
Anderson collective phonon mode appearing as a pole in x- It is evident that the origin 
of this pole lies in the vertex correction, since this is also the pole of T^ 2 \ The pole is 
given by 

<[u; 2 -(v F .p 9 ) 2 ]/> = 0. (68) 
In the limit q — > and u — > 0, / ~ 1 and hence the pole is 

^BA = ~^ V FPl ( 69 ) 

The BA mode restores the continuous symmetry which is broken by BCS ground state. 
It is required to fulfill the Ward identities |3H]- In the low momentum and low energy 
limit (0 < /3 « 1) the dynamic structure factor can be obtained by linearizing the 
denominator of the second term in Eq. ()48j) around the BA mode. By approximating 
/ ~ 1, we then obtain 

S(q, u) = AT(0)73 2 r^-% - <"ba). (70) 
Zujba 

With 73 — > 1, this satisfies the /—sum rule 

Nq 2 

uS(u,q)duj = (71) 

o 2m 

where N is the total number of particles. 

To have higher order (in terms of £g) corrections [5T] to the BA mode, we expand 
the function /(/?) to the fourth order in (3 and obtain the result 

/(/^l + ^T + f^ (72) 
Then the pole is then given by 

resulting in 



{(3 2 f(P)} ^((3 2 + J^} = (73) 



2 _ (y F p q ) 2 



1 _ A ( v fPq\ 
45 V 2A J 



3 - lr 4 ox , (74) 

BA mode is well defined in the low momentum regime, i.e., for £q = VFP q /{2A) « 1. 
For large momentum, it becomes ill defined due to Landau damping. To get the 
dynamical correction to BA mode [33 E2] , the right hand side of Eq. ([74^1 needs to 
be multiplied by a factor [1 — y|a s (iV(0)]. 

Before closing this section, we would like to stress that the polarization-selective 
small angle stimulated light scattering may be useful in exciting BA mode. Because, 
unequal momentum and energy transfer can be accomplished by making 7^ 7^ 7^. 
This will lead to unequal response of the two spin states. In the small momentum and 
energy transfer regime, this will result in large wave-length center-of-mass motion of 
Cooper-pairs and hence superffuid density fluctuation [33] . However, how to detect this 
BA mode of superffuid trapped atoms is presently unknown. 
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Figure 7. Scaled dimensionless dynamic structure factor (DSF) of superfluid trapped 
Fermi g 8jS clS cl function of dimensionless energy transfer (ui/ep) for different cases (see 
the text). For the sake of comparison, DSF for the two cases are scaled differently: 
Si(u), q) and Su(uj, q) are scaled by the factors 1/[27V(0)] and l/iV(0), respectively. 
The momentum transfer is kept fixed at q/kp = 0.2. The scattering length is 
a s = 0.51/c^ 1 for which the BCS gap Eq. H76|) yields the value of the gap at the 
trap center as Ao = O.ObeF- Case-I (dotted) refers to the unpolarized light in the 
absence of magnetic field, case-II refers to the circularly polarized light in the presence 
of magnetic field. Dashed curve (Eq. (57)) is for unpolarized light without vertex 
correction when only BCS-type mean-field is used. 




Figure 8. Same as in Fig. 1 but for q = 0.8kp 



8. Numerical results and discussion 



We now apply the formalism discussed above to harmonically trapped superfluid Fermi 
atoms. For simplicity, we consider an isotropic optical trap characterized by the length 
scale ato = \Jh/ (muho), where u)ho is the trapping frequency. In Thomas-Fermi local 
density approximation (LDA) [53J, the state of the system is governed by 

e F {v)+V ho {r) + U(r)=^ (75) 
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Figure 9. DSF of BCS superfluid (solid) and normal fluid (dashed) are plotted as a 
function of energy transfer for q — OAkp and Aq = 0.05 for the case-I. 



where e^(r) = h 2 kp{r) 2 / [2m) is the local Fermi energy, kp{r) denotes the local 
Fermi momentum which is related to the single-spin local number density by n(r) = 
kp(r) 3 /(6tt 2 ). Here U represents the mean-field interaction energy and \x is the chemical 
potential. At low energy, the mean-field interaction energy depends on the two-body s- 
wave scattering amplitude fo(k) = —a s /(l + ia s k), where a s represents s-wave scattering 
length and k denotes the relative wave number of two colliding particles. In the dilute 
gas limit (\a s \k « 1), U becomes proportional to a s in the form U(r) = ^f^n(r). In 
the unitarity limit \a s \k — > oo, the scattering amplitude / ~ i/k and hence U becomes 
independent of a s . It then follows from a simple dimensional analysis that in this limit, 
U should be proportional to the Fermi energy: U(r) = j3 u ep{v) where (3 U is the constant. 
In this limit, the pairing gap also becomes proportional to the Fermi energy. Based on 
the regularized mean-field approach discussed earlier and LDA, the zero-temperature 
density profiles [S3], momentum distribution |55j and the finite temperature effects [SB] 
of superfluid trapped Fermi atoms have been recently studied. For dilute gas limit, the 
local density distribution of trapped gas may be may be approximated by neglecting the 
interaction term U in Eq. (|T5]). In the BCS limit {kpa s — > 0~), the gap is exponentially 
small and can be expressed by the well known formula 



A B cs - —5- exp 



TV 



(76) 



where ef is the Fermi energy. 

Under LDA, the density profile of a trapped Fermi gas is given by 

n(r) = n(0)(l-r 2 /^ F ) 3 / 2 , (77) 

where n(0) = 1/ (6ir 2 h 3 )[2m[j,/ (1 + P u )} 3 ^ 2 is the density of the atoms at the trap center. 
Here R^ F = 2/i/ (mw^) is the Thomas-Fermi radius of the trapped atomic gas. The 
normalization condition on Eq. (f77|) gives an expression for [i = (1 +/3 n ) 1 / 2 (6iV CT ) 1 / 3 fru;/ lo 
where N a is the total number of atoms in the hyperfine spin a. The Fermi momentum 
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Figure 10. Dimensionless scaled DSF 5/(w,q)/(2iV(0)) (dashed) and Sn((J,q)/N(0) 
Vs. dimensionless energy transfer oj/ep for q = OAkp and a s = 3.21/c^ 1 . For this 
larger scattering length, using the regularized gap Eq. J2fJJ|, the gap is calculated to 
be A = 0.15e_p. The inset shows the corresponding plot in Case-I for leading order 
approximation as given Eq. (57). 



at the trap center k% = [37r 2 n(0)] 1/3 = (1 + f3 u )~ 1/A k F where k F = (48A^ (J ) 1/6 / 'a ho is the 
Fermi momentum of the noninteracting trapped gas. Under LDA, the dynamic structure 
factor is given by 



where 5 r (q, u) is the DSF for Fermi momentum kp{r) evaluated at a position r assuming 
the system is locally uniform. Here Vtf = (4:/3)ttR ff is the Thomas-Fermi volume. 

In Figs. I71ITIH we show DSF (calculated using LDA) of superfluid trapped Fermi gas 
due to single-particle excitations only as a function of energy transfer under different 
physical conditions. Figs. OElare plotted for a s k F < 1 with gap given by BCS gap Eq. 
(I75J) while Fig. [TU] is for a s k F > 1 with gap determined by the regularized gap Eq. (|31j|) 
coupled with the superfluid number Eq. (j3*Tj) . For the sake of better comparison, we have 
scaled the DSF in Case-II (denoted by Su hereafter) by a factor [iy(0)] _1 while DSF in 
Case-I (denoted by Sj hereafter) is scaled by half of this factor, that is, by [2iV(0)] _:L . 
Here N(0) refers to the density of states at the trap center. By comparing Fig. 0and Fig. 
IH] which are plotted for lower and higher momentum transfer, respectively, we infer that 
the vertex correction is most significant in low momentum and energy transfer regime. 
At high momentum and energy transfer regime, mean-field approximation seems to be 
reasonably good. Furthermore, at lower momentum transfer, Sn shows larger deviation 
from Si/2 with both tending to equalize at higher energy transfer. 

In Fig. El we compare DSF of superfluid gas with corresponding DSF for normal 
fluid. The two curves do not show any discernible shift of their peak values apparently 
due to exponentially small gap. However, as u> decreases below the value at which 
the maxima occurs, DSF in superfluid case exhibits reducing gradient in contrast to 
that of normal case of almost steady gradient. This feature may constitute an indication 




(78) 



Light scattering in Cooper-paired Fermi atoms 



22 



of the occurrence of BCS-type superfluidity in trapped Fermi gas. This feature can 
be explained on the basis of inhomogeneous density distribution of trapped gas. For a 
uniform Fermi superfluid, in the single-particle excitation regime, DSF remains zero until 
energy transfer exceeds 2A at which it rises sharply with the increasing energy transfer. 
For a superfluid trapped Fermi gas, owing to the spatial distribution of the gap, DSF 
has a structure below 2A , where A represents the gap at the trap center. As to goes to 
zero, the gradient of S(5, q) vanishes. In the low energy regime {u < 2A ), S(u, q) varies 
with uj nonlinear ly. When uj approaches 2Ao, the gradient changes abruptly implying 
a discontinuity (which may be indiscernible experimentally on practical grounds). This 
behavior can also be explained by considering the boundary condition 2A(x) < u. This 
spatially dependent lower bound on u implies that, when uj is less than 2A , the atoms 
at the central region of the trap can not respond to the light fields via single-particle 
excitations, only those atoms in the peripheral region can do so. 

Fig. El displays DSF for both the Case-I and Case-II for larger scattering length 
a s = 3.21/c^ 1 for which the gap is A = 0.15ei?. For both the Figs. El and EU q is 
fixed at OAkp. In comparison to the Fig. El we notice that the peak of DSF in Fig. fTOl 
exhibits a shift apparently due to the occurrence of relatively larger gap. We further 
notice that the width has been broader with peak value reduced by roughly one order of 
magnitude. This may be attributed to the relatively larger interaction and hence larger 
vertex correction. 

9. conclusion 

In conclusion we have presented a detailed theoretical analysis of the response of Cooper- 
paired Fermi atoms due to off-resonant light scattering at zero temperature. We have 
studied vertex correction which is quite significant at low momentum. By making use of 
the Zeeman shifts between two ground hyperfine spin states and also between the excited 
state hyperfine spin manifolds, we have shown that it is possible to transfer different 
amount of momentum and energy to the two partner atoms of a Cooper pair. Light 
polarization plays an important role in selective single-particle excitations in superfluid 
Fermi atoms. Using circularly polarized light in the presence of a magnetic field, quasi- 
particle excitation can be obtained in one spin component only. We have analyzed 
the dynamic structure factor (DSF) due to single-particle excitations under different 
physical conditions. DSF shows a shift for large gap. In contrast to trapped normal 
Fermi atoms, the gap inhomogeneity of trapped Cooper-paired Fermi atoms leads to 
relatively reduced gradient of DSF below 2A , where A is the gap at the trap center. 
This reducing gradient may constitute a signature of superfluid state. 

Although the focal theme of this paper has been the theoretical analysis of 
stimulated scattering of polarized light by superfluid Fermi atoms under different 
physical conditions, there is some relevance of it in experiments with Fermi atoms. 
The question arises how to detect experimentally spectrum of density fluctuation (or 
DSF) of superfluid two-component Fermi atoms such as 6 Li. Towards this end, we wish 
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to present some speculative and suggestive discussions. We recall that the DSF of BEC 
has been experimentally detected |SH] using stimulated light scattering (or Bragg 
spectroscopy) and the well-established method of time-of-flight measurements. One of 
the major difficulties in evaluating DSF of two-component Fermi atoms from time-of- 
flight measurements might stem from the fact that the initial information about the 
momentum distribution of the atoms may be washed away during expansion due to 
relatively large s-wave collisions of the two hyperfine spin components. This difficulty 
can be circumvented by the method of rapidly reducing the magnetic field (that induces 
Feshbach resonance) just before switching off the trap as done in numerous recent 
experiments [HI EJ 121] • In light scattering experiment, polarization-selective light 
scattering may be useful in suppressing the collisions among the scattered atoms during 
their expansion on being released from the trap. Since the scattered atoms will be in 
a single spin state, there will be diminished probability of collision among those atoms 
(the leading order p-wave collision at low temperature is vanishingly small). This may 
lead to better precision in time-of-flight spin-selective measurements |2SIEI| of scattered 
atoms. Order of magnitude analysis of Ref. [12] suggests that, with large momentum 
transfer, it may be possible to distinguish the scattered atoms in time of flight images. 
To reveal the information about the momentum and density distribution of scattered 
atoms, the time-of-flight images with and without Bragg pulses should be compared. 
Furthermore, Bragg spectroscopy allows one to choose different directions for scattered 
atoms, since the scattering is of predominantly stimulated type. It may be possible to 
scatter atoms in two opposite directions by using three or four beam stimulated light 
scattering configuration as discussed in One can then explore the possibility of 

measuring the correlation of two scattered atoms with opposite momentum by similar 
technique as applied in recent theoretical and experimental [SU] studies. Finally, 
polarization-selective light scattering may be useful in exciting BA mode the detection 
of which poses a challenging experimental problem. 

Appendix 

We here present an alternative derivation of DSF without vertex correction using BCS- 
Bogoliubov mean-field treatment. This DSF coincides with that obtained using leading 
approximation as described in the text. It can be defined by S(q, uj) = J2f\(f I 
Z)t7=T,j. 7<«>-Pct( < i) I 0)| 2 5(cj — €f + e ) where |0) represents the many-body ground state 
with energy eo and the sum runs over all the final states |/) which can be coupled to the 
ground state by the density operator p CT (q) = J2k a t,k+q a °-,k- The DSF for the Case-I, 
that is, for the condition = has been explicitly calculated in Ref. jHZj. We here 
assume ^ 7j_j_, and as an extreme case we consider the case-II: 7^ ~ and — > 1. 
Setting e = 0, we can write 




(A.l) 
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where k' = |k + q| is the wave number of a scattered atom, V is the volume of the 
system and — v\ — (1 — £,k/Ek)/2 is the momentum distribution function. Here 
Ek = \Ji1 + A 2 is the energy of an elementary excitation (Bogoliubov's quasiparticle), 
A k represents the pairing gap and = h 2 k 2 / (2m) — fx. Note that the usual BCS 
coherence factor [SHIS] m(k,k') = Uk'Vk +UkVk> where %i\ = 1 — v%, has changed. This 
is due to the fact that the momentum and energy transfer occurs in either partner of a 
Cooper pair because of polarization-selective light scattering. 

We here give an outline of the method of calculation of the integral in Eq. (jA.l|) . 
For notational simplicity, we denote C — Cfc an d £' = 6fc' — £ + Vfc-Pg! and similarly we 
replace Ek and E^ by E and E', respectively. The integration may be restricted near 
£ = €k — fi ~ 0, since the dominant contribution to the integration comes from /c-values 
near ~ /i. For convenience, we change the variable of integration into E by using the 
relation d£ = EdE/(E 2 - A 2 ) 1 / 2 . Using the identity 

= Hi (A ' 2) 

where E is the solution of the equation E + E' = u, we have 



( E -o(E'+e) 



(A.3) 



lE=E 



4 | 1 + dE'/dE | 
After a lengthy algebra as in Ref. [37], we then obtain 

where /? is defined in Eq. (I28|) . Let z = cos#, where 9 is the angle between k and q. 
Changing the variable of angular integration into x = VkP q z / \/uj 2 — 4 A 2 , one obtains 

2p q v F uj Jo (1 - jx 2 ) 3 / 2 Vl - ^ 

where j = 1 — 4A 2 /cj 2 and x = Min 1, r^z^^yjz • m writing the above equation, we 
have considered the weak-coupling case and hence replaced the chemical potential \i by 
the Fermi energy ep. If 2A < uo < (p q VF) 2 + 4A 2 ) 1//2 , then xq = 1 and the result is 

5(q, 6) = + M 2 iq(3/2, 3/2; 2, j)] (A.6) 

where M = 7rj(l — j)/4. Here -E(j) represents the complete elliptic integral and 
2^1(0, b; c, d) is the hypergeometric function. In the limit A — > 0, ^(j) — > 1 and 
M 2 F\ — > 1 leading to the result 5(5, q) = / '{ApqVp) which is half the dynamic 
structure factor of normal fluid jlH] consisting of noninteracting quasi-particles within 
the specified energy range. The factor half arises due to our initial assumption 7^ ~ 
and 7j_j_ — > 1 which implies that only spin j are scattered and since the number of atoms 
in the two spin components are assumed to be equal, only half of the total number of 
atoms contribute to the scattered flux of atoms. 
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